Stochastic models for a chemostat and long 

time behavior 

P. Collet* Servet Martinez} Sylvie MeleardJ Jaime San Martin^ 

March 8, 2013 

Abstract 

We introduce two stochastic chemostat models consisting in a cou- 
pled population-nutrient process reflecting the interaction between the 
nutrient and the bacterias in the chemostat with finite volume. The 
nutrient concentration evolves continuously but depending on the pop- 
ulation size, while the population size is a birth and death process 
with coefficients depending on time through the nutrient concentra- 
tion. The nutrient is shared by the bacteria and creates a regula- 
tion of the bacterial population size. The latter and the fluctuations 
due to the random births and deaths of individuals make the popu- 
lation go almost surely to extinction. Therefore, we are interested in 
the long time behavior of the bacterial population conditioned to the 
non-extinction. We prove the global existence of the process and its 
almost sure extinction. The existence of quasi-stationary distributions 
is obtained based on a general fixed point argument. Moreover, we 
prove the absolute continuity of the nutrient distribution when con- 
ditioned to a fixed number of individuals and the smoothness of the 
corresponding densities. 
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1 Introduction 



Since some decades and the first work of J. Monod [9] and Novik and Szilar 
[TO] . [TTj . see also p2], biologists have developed procedures which allow to 
maintain a bacterial population at a stationary finite size while at the same 
time, the bacteria have a constant individual growth rate. The procedure is 
based on a chemostat: bacteria live in a growth container of constant volume 
in which liquid is injected continuously. The entering liquid contains a fixed 
concentration of nutrient but no bacteria (fresh liquid). In the container, the 
nutrient is consumed by the bacteria. We assume that the chemostat is well 
stirred, so that the distribution of bacteria and nutrient are spatially uniform. 
Since the container has a finite volume and fresh liquid enters permanently, 
an equal amount of liquid pours out containing both unconsumed nutrients 
and bacteria. 

These chemostats are extremely useful in the study of bacterial population 
dynamics, in particular in the study of the selection of the fastest growing 
species or the fixation of advantageous mutations (see [I], [5J, [B] or [8]). In 
the literature, their study is mainly based on deterministic models where 
both nutrient and bacteria population dynamics are described by a coupled 
deterministic continuous process. Deterministic approximations for the bac- 
terial population's size are justified by a large population assumption. 

In this work, we develop stochastic chemostat models based on a previous 
work of Crump and O' Young [3 J taking into account that the bacteria popu- 
lation may not be large enough so that a deterministic approximation can be 
justified. We introduce a coupled population- nutrient process reflecting the 
interaction between the nutrient and the bacterias in the chemostat. The 
nutrient concentration evolves continuously but depending on the popula- 
tion size, while the population size is a birth and death process with coef- 
ficients depending on time through the nutrient concentration. Moreover, 
the time derivative of the nutrient concentration jumps simultaneously with 
the population size. We thus take into account the random fluctuations of 
this population size due to the individual births and deaths. The bacteria 
need nutrient to reproduce. We will consider two cases. In the first one, the 
bacterial population dies instantaneously if the nutrient is missing. In the 
second case bacteria can survive without nutrient by undergoing some kind 
of hibernation and may wake up once nutrient reappears. The nutrient is 
shared by the bacteria. This creates an indirect competition between bac- 
teria and leads to a regulation of their population size. In our models, the 
fluctuations due to the random births and deaths of individuals and the size 
regulation make the population go almost surely to extinction. Therefore, 
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the long time behavior of the population's size is obvious and the interesting 
questions concern firstly the rate of extinction and secondly the long time 
behavior conditioned to the non-extinction which is captured by the notion 
of quasi-stationary distribution. 

To our knowledge, the models introduced in this paper are the first stochastic 
chemostat models where interaction between bacteria is taken into account 
leading to extinction. The study of quasi- st at ionarity gives nevertheless a 
description of a quasi-stability which can happen in a faster time scale than 
extinction. This work concerns monotype individuals but could be general- 
ized to a multi-type population. 

In Section 2, we describe the two population-nutrient models described above 
and prove in Section 3 their global existence. We also show the extinction 
of the population when time increases. The existence of quasi-stationary 
distributions is obtained in Section 4. Our main theorem is based on a general 
argument proved in [2]. In Section 5, we prove the absolute continuity of the 
nutrient distribution when conditioned to a fixed number of individuals and 
the smoothness of the corresponding densities. 

2 The Stochastic Population-Nutrient Pro- 
cess 

We consider a stochastic discrete population process describing the dynam- 
ics of a bacteria population for which individuals develop and reproduce 
depending on the quantity of nutrient y in the solution. The dynamics of 
the nutrient is related to the consumption of the individuals. We assume 
that the concentration of nutrient in the injected solution (without bacteria) 
is a constant equal to y*. The chemostat has a finite volume equal to one. 
The liquid enters in the chemostat free of bacteria and pours out after being 
well stirred in the container. The pouring out liquid contains bacteria. The 
dilution coefficient of nutrient in the fresh liquid per unit of time is D. Since 
the liquid is well stirred, around N(t)D bacterias will be washed out in the 
pouring out, when N(t) is the size of the bacteria population at time t. Thus 
D is also the rate at which an individual will disappear due to the evacuation 
of liquid. 

We consider coupled processes in which the nutrient concentration evolves 
continuously while the bacteria population size evolves as a time-continuous 
birth and death process with coefficients depending on the nutrient concen- 
tration. We assume that the nutrient is partly consumed during the repro- 
duction of bacteria. 
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We will denote by (Y(t),t > 0) the concentration of nutrient and by (N(t) : 
t > 0) the population size process. The stochastic process Z = (Z(t) := 
(N(t),Y(t)) : t > 0) describes both the population size and the nutrient 
concentration in the chemostat. 

Let us now define the parameters of the model. 

If y is the quantity of nutrient, then the birth and death parameters driving 
the dynamics of the population are as follows. 

• The birth rate per individual is b(y), where the function b : M + — > R + is 
assumed to be an increasing continuous function and such that 6(0) = and 
b(y) > for y > 0. We assume that b is bounded with an upper-bound b^ . 

A usual example of a function b is, for some constant K > 0, 

y 



b(y) = &oc 



K + y 



An extra hypothesis that we will add for some results is that b is differentiable 
and f (0) > 0. 

• The background death rate per individual is d(y), so it is supposed to be a 
function of the concentration of nutrient. The function d : K + — > IR + U {oo} 
is assumed to be continuous non-increasing, strictly positive, and d(0) is the 
unique value that can be infinite. 

• The dilution makes each individual disappear at rate D independently of 
the birth and death events. 

• The per individual rate of consumption of nutrient for reproduction is 
^r, where R denotes the biomass yield. Furthermore, individuals consume 
nutrient during their life and the quantity of nutrient consumed per individual 
will be denoted by rj > 0. 

• We will consider two cases. 

In the first case, individuals need nutrient to survive. Then we will assume 
that their death is instantaneous as soon as nutrient is missing, therefore 
d(0) = oo. 

In the second case, bacteria enter in some kind of hibernation if nutrient is 
missing. That means that c?o(0) can be finite. 
In both cases, we will set 
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Let us now describe the process. In both cases the nutrient concentration 
process Y = (Y(t) : t > 0) evolves according to 

^ = D(y*-Y(t))-b(Y(t))N(t). (1) 

The process Z has the following infinitesimal generator: for (n, y) G N x M. + , 

Cf(n,y) = b(y)nf(n + l,y) + (D + d(y))nf(n-l,y) 
-(b(y) + D + d(y))nf(n,y) 

+ ( D (y*-y)-nb(y)} d y f(n,y) . (2) 

We refer to |TJ for the numerical study of similar models. 

Below, we show that the hypotheses we gave on the coefficients leading the 
process, guarantee that Z is well-defined and takes values in N x [0, y*]. 

Proposition 2.1. The process Z is well defined and takes values inN x M + 
for all positive time t G K+. Moreover, N + x [0, y*\ is an invariant set for 
the process Z, so Y(0) G [0,y*] implies Y(t) G [0,y*] for all t G R + . 

Proof. The process iV = (N(t) : t > 0) has a pathwise representation driven 
by a Point Poisson measure Af(dQ, ds) defined on R + x M + : 



N(t) 



N + J l e < b{Ys) N{d9,ds)- J l b (Y 3 )<9<b(Y 8 )+D+d(Y 3 )N{d9,ds); 



E+x(Q,t] R+x(0,t] 

F(0 = y(0)+jTr£>(y*-y(s))-6(Y(s))iV^d s . (3) 

It is obvious that the process iV is stochastically upper-bounded by a birth 
process with individual birth rate boo. This latter does not explode, so does 
N. 

Let us now study the nutrient concentration Y = (Y(t) : t > 0). Note that 
the assumptions on the parameters and Equation ([T]) ensure that Y(t) > 0. 
Indeed, the derivative of Y(t) at y = cannot be negative. 

Let us show that N x [0, y*\ is invariant. Take Y(0) G [0,y*]. A standard 
comparison theorem yields Y(t) < v(t) where v'(t) — D(y* — v(t)) ; v(0) = 
Y(0). But in that case, v(t) = y* — (y* — Y(0))e~ Dt < y* then it remains in 
[0, y*} forever proving the invariance. 

□ 

For the initial condition N(0) G N* and Y(0) G R+ \ [0,y*\ we have that at 
time t = (Y(0) — y*)/b(y*) the process has already attained the invariant set 
[0,?/*] x M + or became extinct. 
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3 Study of Extinction 



We are now interested in studying the extinction of the population or the 
complete consumption of the nutrient or other specific states of the population- 
nutrient process and the associated hitting times. 

Let B(N x R + ) be the class of Borel sets of N x R + and B{R+) the class of 
Borel sets of R+. For D e B(N x E+) we put by T D = inf {t > : Z(t) e D} 
the hitting time of D by the process, with the usual convention oo = inf 0. 
The hitting time of the boundaries will be denoted in such a way that the 
reference to process N will be avoided, the contrary for Y where we will 
explicit it. More precisely, 

To = T {o}xM +=inf | t > . AT(t)=0}, T y=0 = T Nx W=inf{t>0 : Y(t) = 0}. 

(4) 

We will also denote by T< m = T^°'-' m ^ xR + = inf{t > : N(t) < m}. Analo- 
gously for T <m . 

Note that the set {0} x M + is an absorbing set, that is N(t) = for all t > T . 

After T the nutrient Y(t) is absorbed linearly at y*, in fact: Y(t) = y* 
for t > T + (y* - Y(T ))/Dy* and Y(t) = Y(T ) + D(t - T ) for t G 
[T ,T Q + {y*-Y(T ))/D]. 

Let us see what happens after Ty =0 . In the case rj > Dy* the nutrient 
Y(t) remains at up to the extinction of the population, so Y(t) = when 
Ty =0 < t <T . The evolution of the nutrient after T was already described. 
Now assume nrj > Dy* for some n > 1, and denote hq the minimal of 
these values, so (no — 1) rj < Dy* < n rj. In this case Y(t) = when 
T Y =o < t < T <no . In the case where d (0) = oo, we get 

Ty=o < T =^ Ty =0 = T 
because all individuals die instantaneously. 

Let us firstly study the stationary nutrient concentration states at fixed pop- 
ulation size. 

Lemma 3.1. Consider the equation 

G n (y) = with G n (y) :=D(y*-y)-n b(y) . (5) 

Then, 

{i) If r] = 0, then for any neN, Equation ([3]) has a unique simple root y n , 
which belongs to [0,y*]. In addition, the sequence (y n : n 6 N*) of the roots 
decreases to 0. 
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(ii) If r] > then Equation fl5J) has no root for n > and admits a simple 
root y„ for n < . 

Proof. Let us fix n. Obviously for n — 0, Equation ([5]) has the trivial root 
y = y*, so we restrict n to be in N*. By assumption, the function G n is 
strictly decreasing so for each n E W there exists at most one root. Note 
that for all such n we have G n (y*) < and so there is no root to (J5|) in the 
set [y*, oo). 

Assume rj = 0. For all n we have G n (0) > then there exists a unique root 
which is denoted y n , so it satisfies G n (y n ) = 0. On the other hand, we have 

G n+X (y n ) = D(y* - y n ) - (n + 1) b(y n ) = - b(y n ) . 

Then, G n+1 (0) > 0,G n+1 {y n ) < 0. We deduce < y n+1 < y n < y*. Let 
l/oo = lim y n . By continuity 

n— >oo 

- Z/oo) -nb^yoo) = lim - j/ n ) - nb(y n ) = 0. 

n— >oo 

Then, necessarily 6(j/oo) = and so y^ = 0, and b(y n ) ~ g ^ R as n — >■ oo. 

Assume 77 > 0. Then for all n > Dy* /r] we have G n (0) < and so (jSJ) has no 
solution. Hence the same argument as before gives the existence of a finite 
set of roots (y n : 1 < n < [Dy*/r]\) decreasing with n, where [Dy*/r]\ the 
biggest integer that is smaller or equal to Dy* /rj. 

Note that when n = Dy* /rj then y m = and if Dy* < rj, there is no 
root. □ 



From Lemma 13.11 we know that the set 

S = {y E R+ : 3nE W, D(y* - y) - b(y)n = 0}. (6) 

is a countable set included in [0, y*). If rj = 0, it is infinite and accumulates 
at and if 77 > 0, it is finite. 

In the sequel, when we refer to y n , we will assume implicitly that it exists, 
namely we are in the case rj = or rj > but n e {1, .. \_Dy*/r)\}. 

Corollary 3.2. (i) The set N x [0, y%] is invariant for the process 

(Z(t) :t<T ) up to extinction, that is if Z(0) = (N(0),Y(0)) 6Nx [0, yi ] 

then Z(t) = (N(t), Y(t)) G N x [0, y x ] for allt<T . 

(ii) The set N x [0, y n ] is invariant for the process Z T<n = (Z(t) : t < T <n ). 
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Proof. Let us show the first part. We have Z(t) = (N(t), Y(t)) eW x [0, y x ] 
for all t < T because dY(t)/dt < when y > y\ and n > 1, and so if the 
trajectory arrives to y% the variable Y(t) immediately decreases. For n = 
it is evident, because we stop the process at T , and so Y(T ) = Y(T Q ~) but 
Y(T ) < yi since N(T ) = 1. 

Part (ii) is shown in a similar way as (i). □ 



Let us state a useful lemma. 

Lemma 3.3. For any 6 N*, there exists t > swc/i that 

inf p (n>tf) (r <to) >o. 

ve[Q,i/*],l<»<no 

Proof. This follows at once from the fact that the population process is 
stochastically dominated by a birth and death process with birth rate b(y*) 
and death rate D + d(y*). □ 

Theorem 3.4. We have extinction of the population almost surely, namely 
for any y G [0, y*] and any integer n, 

P(„ )2/) (Tq < oo) = 1 . 

Proof. Remark that an obvious comparison theorem as used in the previous 
proof cannot be applied. Indeed the birth and death rates b(y*) and D + d(y*) 
could correspond to a supercritical case. The effect of the chemostat through 
the nutrient is a regulation of the population. 

We will exhibit an integer n such that the population process will spend an 
infinite amount of time below n + 1. 

Let y > be such that b(y) < D + d(y). Note that by monotonicity, for all 
y E [0, y] we have b(y) < D + d(y). 

Let us define n as an integer such that y no < y if some exists or equal to 

[Dy*/b(0)] + 1 otherwise. 

Let r be the random time defined by 

sup N(t) < n . 

t>T 

Then, it follows from Lemma 13.31 and the Markov property that if P( n ,j/) ( r < 
oo) > 0, then 

^(n,y) {To < OO | T < OO) = 1 . 
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Assume now P( n)2/ )(r = oo) = 1. Let r' be the random time denned by 

inf N(t) > n . 

t>r' 

Assume P( n ,y) (t' < oo) > 0. Let yo(t) be the solution of the differential 
equation 

^ = D(y* - y (t)) - (n Q + l)b(y (t)) , 
with initial condition yo(0) = y*. Let t\ > be the finite solution of 

Voih) = y ■ 

The time t\ is finite from the choice of n . It is easy to verify that for any 
integer valued measurable function n(t) > n + 1, the solution y{t) of the 
equation 

f f = D(y* - y(t)) - n(t)b(y(t)) , 
with initial condition y(0) < y*, satisfies 

y(t) < y 

for any t > t\. On the set { r' < cxo}, the process (N(t),t > r' + ti) 
is dominated by a linear birth and death process with birth rate b(y) and 
death rate D + d^y) (from the monotonicity of the functions). This birth and 
death chain attains no almost surely in finite time since b(y) < D + do(y), see 
[T3] . Hence on { r' < oo }, the process (N(t),t > r' + t\) should also attains 
n in finite time. This contradicts our assumption P( n ,y) (t' < oo) > 0. 
It remains to consider the case r = r' = oo almost surely. In this case there 
exist two infinite sequences of random times 

Ti < S 1 < T 2 < S 2 < ■ ■ ■ 

such that 

N(t) < n for t G [Tj, S d ) , N(t) > n for t E [Sj, T j+1 ) . 

Since we visit the set {iV < no} infinitely many times, and at each visit we 
have a uniformly positive probability of extinction, it follows by the Markov 
property and the Borel Cantelli Lemma that 

P(n,y) {To < OO | T = T = Oo) = 1 . 

□ 

One of our main objectives of this work is to study the processes up to 
the moment the population is extinct Z T ° = (Z(t) : t < To), or before the 
moment of extinction Z T o = (Z(t) : t < T ). All the statements related to 
quasi-stationary distributions depend on Z T o . 
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4 Existence of Quasi-Stationary Distributions 

A quasi-stationary distribution v (with respect to the absorbing time To) is 
a probability measure defined on N* x M + that verifies 

K(Z(t) E D\Tq > t) = v{D) , V£> G B{W x R+) . (7) 

It is known that starting from a quasi- stationary distribution, the time of 
absorption is exponential, that is P„(T > t) — e~ xt where A = X(u) > 0. Let 
Ait, = {/ : N* x M + — )■ R bounded and measurable}. Equation ([7]) can be 
written as, 

POO 

3 A > : K(f(Z(t)), T >t) = e- xt / f(y)u(dy) Vt > 0, / G -M b . (8) 

./o 

Denote by 

Kn = z/({ri} x R+) = P v (JV(0) = n) , 
and the probability measure conditioned to have n individuals by 

u n (B) = is({n} xB|{n}x K+) = P l/ (y(0) G 5 | N(0) = n) MB G . 

Then = £ K n z/„,(rJn{n} xM + ), z/ n (R + ) = 1 VneN* and k„ = 1. 

rieN* rt£N* 

In order that the probability measure v is a quasi- stationary distribution, it 
must satisfy the infinitesimal condition deduced from (18~1), which is given by, 



;>oo 

3 A > : K n / tfe/ n (j/) 



£f(n,y) - A/(n,y) 



V/ G M 6 . 



In terms of the adjoint operator £t of £ defined in (jSJ) the equation to be a 
quasi-stationary distribution is, 

3A>0: C)v = -\v. (9) 



Theorem 4.1. Assume there exists < o < 1 suc/i i/iai limsup?/ "(i(y) < oo. 

Then there exists a quasi- stationary distribution. Moreover there exists at 
least a quasi- stationary distribution such that Y is supported in [0, y*]. 

Proof. It suffices to show the existence of a quasi- stationary distribution. In 
fact, the last part of the statement follows when the existence proof is applied 
to the process Z taking values in the invariant set N x [0, ?/*]. 
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The idea for showing the existence is to use the abstract Theorem 4.2 proved 
in [2J. We assume y G [0,?/*]. We define a function 



Ln>l 



Then a simple computation (it is the same computation as for a birth and 
death process since (fx does not depend on y) leads to 



&pi(y,ri) 



if n > 1 , 
-D - d(y) if n = 1 . 



Therefore if <i is bounded above by a constant d$, then 

e-^+^Vi < e*V < ^ . 
If c? is not bounded, we will prove the lower bound 

inf F nom (N(l) > 1) > Q . 

«0>1,S/oGJO,j/*J 

for a constant Q > 0. By the Markov property it suffices to prove that there 
exists Q > such that 

inf P liJro (JV(t) =1, V0<t<l)>Q. 

2/o£[0,j/*] 

Since iV = 1 on the whole time interval [0,1], y(t) satisfies the differential 
equation 

^ = D( y *-y)-b(y), (10) 

with initial condition y . Since there neither birth, nor death on the time 
interval [0, 1] we get 



F 1>yo (N(t) = 1, V0 < t < 1) = e 



- Jo {d{y(t))+b(y(t)) dt 



Since b(y) is bounded uniformly in y, the above quantity does not vanish if 

d(y(t)) dt < oo . 



It is easy to show that there is a constant c > such that for any y* > y$ > 0, 
we have y{t) > ct for any t G [0, 1]. Since d{y) < 0(l)y~ a for y small, we get 



d{y(t)) dt < 0{l) / y{t)- a dt < oo 
11 



It now follows immediately that 

Q ¥>i < e c (pi < tpx . 

The second function is 

^ 2 (2/,n) = l„>i e a ^ n , 

with a(y) = ay + a , a > and a > 0. A simple computation for n > 1 (it 
is the same computation as for a birth and death process since (f 2 does not 
depend on y) leads to, 

C<f 2 (y,n)=ne a ^ n E(n,y) , 

where 

~(n,y) = b(y) (e a ^-l- ^ + (D+d(y)) (e~ a ^ -l)+Da(y*-y)-a V l y>0 . 

Let us show that there exists A > and Nq such that for all n > Nq, 

sup H(n, y) < —A. 
2/e[o,y] 

Define = (D + (e~^-l) +Da{y* - y) - arjl y>0 . We choose a 
and a such that D (e~ a ° - 1) +Day* < 0. Then C(0) < 0. It follows that 
\/y E [0,|/*], C(y) < -^4 for some A > 0. 

Consider N such that e a(s/ * ) - 1 - ^ < 0. We still have for n > N that 
e a(y*) _ i _ ^ < o. 

Then for n > N , we get S(n, y) < —A. Therefore for any C > there exists 
iV(C) such that for any n > N(C), 

C(p 2 (y,n) < -C<p 2 (y,n) . 

Therefore, for any C > there exists r(C) > (finite) such that 

jTip 2 <-C tp 2 + T(C) ip x . 

(the estimate for n — 1 is by direct computation taking T(C) adequately 
large enough). Hence, 

e"V 2 < e' tc ^2 + • 

In the case where d is bounded above by d* < oo, we now choose C > D + d* 
and t = 1 and apply Theorem 4.2 in [2J. 

In the case d unbounded, we choose C > — logQ and apply Theorem 4.2 in 

□ 
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5 Properties of the Quasi-Stationary Distri- 
butions 

Proposition 5.1. Any quasi-limiting distribution has support in N* x [0, ?/*]. 

Proof. We recall that a quasi-limiting distribution v is a probability measure 
on N* x R + such that for some initial point (n 0} y ) efx W + 

u(A) = lim F (no . yo) (Z(t) G A\T > t) , VAe B(N* x R + ). 

t— >-oo 

If |/o G [0,?/*], the assertion follows from Proposition 12.11 Let us now assume 
that yo > y*. We introduce the function t — > v(t) defined by 

-(t) = D(y*-v)-b(y*) ;v(0) = y . 

Let T (y ) be defined by v(r (y )) = y*. For t < T A r (y ), we have 
Y(t) < v{t). It follows from Proposition O that P (no , W)) (Y(t) > y*\T > 
t) — , Vt > To (3/0)5 which concludes the proof. □ 



Theorem 5.2. For a// n <E W and any quasi- stationary distribution, the 
probability measure v n is absolutely continuous with respect to the Lebesgue 
measure, with C°° -density on the set R + \ {0,y n }. 

The proof of this theorem is obtained from the following lemmas. Recall that 
the set S has been defined in (J6]). 

Lemma 5.3. For all n G N* ; the measure v n satisfies 

dv n = c%5 + ( c i 5 vj + U ^V) d v) ( n ) 

with c™ > for j G N and u n is the density of the absolutely continuous part 
of v n (so it is a non-negative integrable function) and it is a C°° function 
outside S U {0}. 

Proof. The measures v n satisfy in the sense of distributions the differential 
equations 

d y ((D(y*-y)-nb{yj) v n ) 

= b(y)nv n -i+(D+d(y))nv n+1 -(b(y)+D+d(y))nv n +\v n . (12) 

Since the right hand side is a measure, we conclude by a recursive argument 
that the measures v n have a C°° density on (S U {0}) c . This also shows 
that these measures have no singular part, and the Lebesgue decomposition 
theorem gives relation (fTTjl . □ 
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Lemma 5.4. Let I be an open interval included in (<SU{0}) c . If there exists 
neW such that u n (I) = 0, then = for all j G N*. 

Proof. Assume that there is an interval I not intersecting S U {0} such that 
for some integer n G N*, v n {T) = vanishes on I. From f[T2j) we deduce 
that u n+ i and v n -i also vanish at /. Therefore, we conclude recursively that 
= for all j G N*. □ 



Lemma 5.5. The probability measure v n is absolutely continuous on M + \ 
{0,y n } and its density is bounded on any compact set contained in M + \ 
{0,2/4- 

Proof. Let us show that if j ^ n, v n cannot have a Dirac mass in yj. We 
will do it by contradiction, so assume it does. Let / be a C°° function with 
compact support containing yj and such that its support does not contain 
any other point of S U {0} except yj. By using formula (TTTT) we get, 

-f'( yj ) (p{y*- yj )-nb(yj) - Jf(y) (p(y* -y) - nb{y)) ^ n {y)dy 
= J f(y) (%) n dv n ^i{y) + (D + d(y)) n du n+1 (y)) 
+ / f(y) (A - (b(y) + D + d{y)) n) dv n (y) . (13) 



It is not difficult to construct a sequence of functions (f q : q E N*) contained 
in C°° with support in a fixed small enough neighborhood of yj and such 
that, 




This leads to a contradiction, when we take / = f g in (1T31) and make q tend 
to infinity. 

Now, using again equation (|T2|) . we deduce easily the boundedness of the 
density of v n outside a neighborhood of {y n , 0}. □ 

Lemma 5.6. The probability measure v n cannot have a Dirac mass in y n . 

Proof. Assume it does. Then in a neighborhood of y n , by using ffTT]) since 
Cj = for j = 0, n and by writing c n = c™, we can write 

dv n = Cg5 + c n 5 yn + u n {y) dy . 
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Let / be a function C°° with compact support containing y n but that does 
not contain 0. We have 



- J f'(y) (D(y*-y)-nb(y)) u n (y)dy 

f(y) {Ky) n du >n-i{y) + (D + d(y))ndv n+1 (y)) 



+ J f(y) (A - (6(1/) + D + %)) n) u n (j/) A/ 

+Cn /(i/ n ) (A - (b(y n ) + D + d{y n )) n) . 

We now construct a sequence (f q '■ q € N*) of C°° functions with support in 
a fixed small enough neighborhood of y n such that for some constant C, 



f q (Vn) = 1 Vg ; sup |/ g (y)| = 1 ; lim |/ g (j/)| = ^ y n 

q— >oo 



sup |y - y n \\f' q (y)\ < C Vg ; lim (y - y n )f'{y) = Vy ^ y n . 



and 



Such a sequence can be easily constructed. 

Recall that v n -i an d zVi-i are absolutely continuous with C°° density near 
y n . We conclude that if c n ^ 

A = (b(y n ) + D + d(y n ))n . 

In the case d is constant, it is known that the eigenvalue A satisfies A = 
n n (D + d) < D + d, and so it is strictly less than (b(y n ) + D + d(y n )) n and 
we obtain a contradiction. 

If d is not constant, the proof of the contradiction is more intricate. 
Coming back to Equation fTT2|) for u n in a neighborhood of y n but outside 
that point we get 

d y (G n (y) u n (y)) = f n (y) G n (y) u n (y) + h n (y) 

with (see Lemma [3. II) 

G n (y) =D(y*-y)-n b(y) = (3 n (y - y n ) + O ({y - y n ) 2 ) , 

with (3 n < 0, and 

,/x %«) ~ b(y) + d(y n ) - rf(y) 

/ " (!/) = " gS) = ' 
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and 

K{y) = b(y) nu n-i(y) + (D + d(y))nu n+1 (y). 

and recall that h n (y) is C°° near y n . The only solution which is integrable 
near y n is given by 

u n (y) = ^yexp f n (s)ds^ J\-fLf^ dw h n (s)ds . 

If h n > on a subset of positive measure of a small neighborhood of y n , we 
have u n < which is a contradiction. 

Therefore, /i n must vanish on both sides of y n . By the above result on the 
support, we conclude that u n vanishes in a neighborhood of y n as well as all 
the Uj with j ^ n (see Lemma l5~il) . In particular, if we consider the equation 
for z/ n+ i in this neighborhood (see ffl2l) ). we get 

= b(y n ) c n 5 yn 

which contradicts c n ^ 0. □ 

Theorem 5.7. On (0, yi), the density of v n satisfies u n > except perhaps 
in y n . 

The proof uses two lemmas. 

Lemma 5.8. If v n +i or v n-\ has a support dense in (0, 2/1) then u n can be 
only in y n . 

Proof. The function G n has a simple zero in y n . Assume that for z G (0, yi), 
z ^ y u u n (z) = (and G n (z) ^ 0). 
Computation as in the proof of Lemma 15.61 gives 

u n (y) = ^y ex P (J* fn ^ ds ) J*e- f:fnMdw h n (s)ds . 

The conclusion follows. □ 

Let us show that the process is irreducible up to extinction and before yi. 

Lemma 5.9. Starting from any initial condition on N* x (0,yi) the process 
Z has dense support on N x (0,yi), that is, 

Vm G N, V?/ G (0,yx), V7 > 0, 3t(y') > such that Vt > t(y'), 
V(yo,no)GN* x (0,yi) : F {n(hyo) (N(t) = m, Y (t) G (y'- 7 , y'+ 7 )) > 0. 
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Proof. In the proof we will assume yo < y' (the case yo G (y', yi) is shown 
similarly). Let > be smaller than min(?/o,2/' — 2/o)/2 and let t > be 
fixed. Then, there is e = e(no, (3) such that the following event has a strictly 
positive probability: 

— On the interval time [0, e] there are exactly n — 1 deaths and there is no 
other jump of N, and so N(t) decreases from n to N(e) = 1; 

— Y(e) belongs to (Y(0) - f3, Y(0) + 

— On the interval of time [e, e + t\ there is no jump of iV (no birth and no 
death) ; 

— On the interval of time [e + t, e + t + 0] there are exactly m — 1 births and 
no other jump when m > 1, there is no jump if m = 1, or there is a unique 
death and no other jump when m = 0; 



is bounded below by D(y* - y') - b{y'). Take t = (y' - Y(0) + P)/{D{y* - 
y') — b(y')), let us see that the number t(y') — e + 1' makes the job. In fact, 
we have ensured that in a time smaller or equal to t(y') we have attained 
{m} x [y' — 7, y' + 7]. For any time bigger than t(y') it suffices to modify 
slightly the above argument and allow a sequence of jumps up to the moment 
that Y(t) has negative derivate and in this way we can postpone the time 
of attaining the set {m} x \y' — 7, y' + 7] from t(y') to a prescribed time 



Proposition 5.10. For alln G N* ; the probability measure v n has a support 
dense in (0, y\). 

Proof. Denote by n ® v n the probability measure defined on N x M + by 
n <8> v n {{fri'} X B) — 5 n (m)v n (B) for all B G S(M+). From Lemma 15.91 and 
the quasi-stationary distribution property ([8]) we have 



Wn® Vn {N(0) = m, y (0) G (y-7, y'+j) = e- xt ¥ n ^ n (N(t) = m, Y(t) G (y'- 7 , y'+j) 




t > t(y'). 



□ 



= e 




P(n, yo )(^W=^ 5 ^We(l/'-7,l/ , + 7))^n(yo) > 0. 



Then, the result follows. 



□ 
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6 Bound on the asymptotic survival rate 

Theorem 6.1. Let A be the exponential extinction rate associated with a 
quasi-stationary distribution v. Then 

A < mf(n(b(y n ) + D + d(y n )). (14) 

n 

Proof. In Section 3 in [2J, it was shown that this extinction rate is an eigen- 
value of the dual problem associated with the probability measure 

C)v = -\v. 

For simplicity, we will prove inequality (I14p in case of n — 1. The general 
case proved in a similar way is left to the reader. 
We introduce the notation 

G(y) = D(y* - y) - b(y) 

H(y)=(D + d(y)) 
The quasi-stationary distribution equation for n = 1 is given by 

Let a e]0, We have for y e]0,yi[ 

4- (G(y)u 1 (y)eS"^ +H ^- x y G ^ ds ) = e^^ +H ^- x ^ ds 2 H(y)u 2 (y) , 



dy 

and integrating between a and y we get 

(G(y) ui(y) eXrW-J+HW-Aj/ow* _ G(a ) Ml(a) 



Therefore 

G(fl)Mi(fl) /«( 6 ( s ) +J? ( s )_A)/G( S )d S 

uiW_ G(j/) 

Let us study more carefully the quantity e~-^( b ( s )+ H ( s )~ A )/ G '( ;s ) d;s . We recall 
that G only vanishes at y\ and that G is decreasing since b is increasing. By 
a simple computation one gets 

m + Jf (.) - A)/G(.) = " ( f' ) + ^,''- A + 0(1) , 
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and thus 

j\b{s) + H(s) - \)/G(s)ds = h{Vl) + G H ^ ] - X log \y - Vl \ + 0{l) . 

Finally, since G(yi) = 0, and G'(yi) < 0, we get for y <y\ 

G(y) = {y- y x ) G\ Vl ) +0((y- yi f) = \y- Vl \ (\G'{ Vl )\ + 0(\y - Vl \)) . 
Therefore 

G(fl)«l(fl) , l -(b(y 1 )+H(y 1 )-X)/G'(y 1 ) 0(1) 

> Gja^a) _ , {b{yi)+H(yi) _ X y G , {yi) g0(1) ^ 
ll/- 2/i I \ G '(yi | 

As we have U\{a) > by Theorem 15. 7\ the integrability of U\ on [0, yi] implies 

6(yi) + g( yi )-A 

G"(m) < ' 

we finally get 

A < b{ yi ) + . 

□ 
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